POINCARE GAUGE THEORY WITH COUPLED EVEN AND ODD PARITY 
DYNAMIC SPIN-0 MODES: DYNAMICAL EQUATIONS FOR ISOTROPIC 
BIANCHI COSMOLOGIES 



Fei-Hung Ho 

Department of Physics, National Central University 
Chungli, 320, Taiwan 

93242010® cc.ncu.edu.tw 

James M. Nester 

Department of Physics, National Central University 
Chungli, 320, Taiwan 
Graduate Institute of Astronomy, National Central University 
Chungli 320, Taiwan 
Center for Mathematics and Theoretical Physics, National Central University 



We are investigating the dynamics of a new Poincare gauge tlieory of gravity model, wliicli lias cross 
coupling between the spin-O"*" and spin-0~ modes. To this end we here consider a very appropriate 
situation — homogeneous-isotropic cosmologies — which is relatively simple, and yet all the modes 
have non-trivial dynamics which reveals physically interesting and possibly observable results. More 
specifically we consider manifestly isotropic Bianchi class A cosmologies; for this case we find an 
effective Lagrangian and Hamiltonian for the dynamical system. The Lagrange equations for these 
models lead to a set of first order equations that are compatible with those found for the FLRW models 
and provide a foundation for further investigations. Typical numerical evolution of these equations 
shows the expected effects of the cross parity coupling. 
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1. Introduction 

All the known fundamental physical interactions can be formulated in a common frame- 
work: as local gauge theories. However the standard theory of gravity, Einstein's general 
relativity (GR), based on the spacetime metric, is a rather unnatural gauge theory. Physi- 
cally (and geometrically) it is reasonable to consider gravity as a gauge theory of the local 
Poincare symmetry of Minkowski spacetime. A theory of gravity based on local spacetime 
geometry gauge symmetry, the q uadratic Poin care gauge theory of gravity (PG, aka PGT) 



search for good dynamical propagating modes led to focusing on the two scalar modes. 

There is no known fundamental reason why the gravitational coupling should respect 
parity. With this in mind, the general quadratic PG theory has recently seen renewed inter- 
est in including all possible odd parity couplings. The appropriate cross parity pseudoscalar 
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was worked out some time 




sketch this theory, and how the 
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coupling constants have been been incorporated, in particular, into the special case of the 
dynamically favored two scalar mode model to give the BHN model," which is the most 
general PG model that we expect to have problem free dynamics. Here we show how a sim- 
ple effective Lagrangian can reveal in a cosmological context the dynamics of this extended 
model. 

We are especially interested in investigating the dynamics of the PG BHN model. This 
can be expected to be very clearly revealed in purely time dependent solutions, hence we 
considered homogeneous cosmologies. The two dynamical connection modes that we wish 
to study carry spin 0+ and spin 0^ (and are thus referred to as scalar modes or more 
specifically as the scalar and pseudoscalar mode). Consequently in a homogeneous sit- 
uation they cannot pick out any spatial direction, and thus they have no interaction with 
spatial anisotropy, so for a study of their dynamics it is sufficient (and most simple) to con- 
fine our attention to isotropic models. Here, we present an extension of the results briefly 
reported in Ref. In particular for the PG BHN model we find an effective Lagrangian 
and Hamiltonian as well as a system of first order dynamical equations for Bianchi class 
A isotropic homogeneous cosmological models and present some sample evolution which 
shows the effect of the cross parity coupling. 



2. The Poincare gauge theory 

In the Poincare gauge theory of gravity,[lEl3]4l5]6l2] {jjg ^-^q ^gj^ of local gauge potentials 
are, for "translations", the orthonormal co-frame '0°- = e"ida;% where the metric is g = 
—1?° (8) 1?" + 5ab'&'^ ® and, for "rotations", the metric-compatible (Lorentz Lie-algebra 
valued) connection 1 -forms T"^ = T^^^^idx"^. The associated field strengths are the torsion 
and curvature 2-forms 

T" :=dz9"+r";3 Ai?'' = ir^^^Tj^Az?", (1) 

which satisfy the respective Bianchi identities: 

DT°' = p A , DR^p-O. (3) 

Turning from kinematics to dynamics, the PG Lagrangian density is generally taken 
to have the standard quadratic Yang-Mills form, which leads to quasi-linear second order 
equations for the gauge potentials. Qualitatively, 

^[d, r] [A + curvature + torsion^] + gi^^curvature^, (4) 

where A is the cosmological constant, k = SttG/c^, and g^^ has the dimensions of ac- 
tion. The field equations, including source terms, obtained by variation w.rt. the two gauge 
potentials have the respective general foms 

A + curvature + D torsion + torsion^ + curvature^ ~ energy-momentum density, (5) 

torsion + D curvature ^ spin density. (6) 
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From these two equations, with the aid of the Bianchi identities (O, one can obtain, respec- 
tively, the conservation of source energy-momentum and angular momentum statements. 

Earlier investigations generally considered models with even parity terms; the mod- 
els had 10 dimensionless scalar coupling constants. The recent BHN investigation" sys- 
tematically considered all the possible odd parity Lagrangian terms, introducing 7 new 
pseudoscalar coupling constants. Not all of these coupling constants are physically in- 
dependent, since there are 3 topological invariants: the (odd parity) Nieh-YarU^ identity 
d(t9" A r„ ) = A Ta + Raf} A ■d°'^ , the (even parity) Euler 4-form R°'P A R''^T]a,3js, and 
the (odd parity) Chem-Symons 4-form i?"^ A R^a- For detailed discussions of the BHN 
Lagrangian and the topological boundary terms see Ref . |8]and the new work Ref . [TT] 

Early PG investigations (especially Refs. [3] fTSl l of the linearized theory identified six 
possible dynamic connection modes; they carry spin-2^, spin-1*, spin-0*. A good dy- 
namic mode should transport positive energy and should not propagate outside the for- 
ward null cone. The linearized investigations found that at most three modes can be si- 
multaneously dynamic; all the acceptable cases were tabulated; many combinations of 
three modes are satisfactory to linear order. Complementing this, the Hamiltonian anal- 
ysis reveale d the related constraints.SSI Then detailed investigations of the Hamiltonian and 
propagatior 1 4 | 1 5 1 1 6 |T7] concluded that effects due to nonlinearities in the constraints could 



be expected to render all of these cases physically unacceptable except for the two "scalar 
modes", carrying spin-0+ and spin-0^. 

In order to further investigate the dynamical possibilities of these PG scalar modes, 
Friedmann-Lemaitre-Robinson-Walker (FLRW) cosmological models were considered. 
Using a A; = model it was found that the 0+ mode naturally couple s to th e accelera- 
tion of the universe and could account for the present day observations] ^ ^ 1 ^ 3 this model 
was then extended to include the mode.^ 

After developing the general odd parity PG theory, in BHN° the two scalar torsion 
mode PG Lagrangian was extended to include the appropriate pseudoscalar coupling con- 
stants that provide cross parity coupling, such terms are often referred to as "parity violat- 
ing" terms. 

o 

The BHN Lagrangian" has the specific form 



1 ^ 

-2A + aoR + boX--Y, ^n^"^^ + icr2V^A'' 



n=l 



1 

'2~g 



12 ^ " 12 ^ + 12^^ 



(7) 



where R is the scalar curvature and X is the pseudoscalar curvature (specifically X/ 6 = 
i?[oi23] is the magnitude of the one independent component of the totally antisymmet- 
ric curvature), and := r"„^ = (^y^^a^, ■= hi^.^^^^cp = ^eM-^^'^^'^T'^a/S 
are the torsion trace and axial vectors. The parameters ag, ai, a-z, W3 and we are scalars, 
whereas 69, 172 and are pseudoscalars. For an extensive discussion of the mathematics 
and physics of the PG theory and this model as well as further references see and 
the new work, Ref. [H] 
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In BHN the general field equations were worked out, and then specialized to find the 
most general 2-scalar mode PG FLRW cosmological model. Here we will take an alterna- 
tive approach and consider manifestly isotropic Bianchi models. 

3. The PGT scalar mode Bianchi I and IX cosmological model 

PG cosmological investigations have a long history. For earUer PG cosmological investiga- 
tions see Minkevich and coworkers, e.g., Refs. i2ni22ll23ll24ll25] and Goenner & Miiller- 
HoissenpS for recent work see Refs. |8]|T8]|T9]|27]|28]|29]|30]|20] 

For the usual FLRW models, although they are actually homogeneous and isotropic, 
the representation is not manifestly so. Indeed they are merely manifestly isotropic-about- 
a-chosen-point. In contrast, the representation used here for the isotropic Bianchi I and IX 
models has the virtue of being manifestly homogeneous and manifestly isotropic. (Indeed 
these are the only two Bianchi models which admit such a representation.) 

For homogeneous, isotropic Bianchi type I and IX (respectively equivalent to FLRW 
k = Q and fc = +1) cosmological models the isotropic orthonormal coframe has the form 

1?° := dt, -d" := acr'^, (8) 

where a = a{t) is the scale factor and cr^ depends on the (not needed here) spatial coordi- 
nates in such a way that 

da' =CeVf^'^^^ (9) 

where C = for Bianchi I and C = 1 for Bianchi IX, thus = fc, the sign of the FLRW 
Riemannian spatial curvature. 

Remark: although a few other Bianchi models can be isotropic, e.g., Bianchi V, the 
representations themselves are not manifestly isotropic, which is an important property in 
our analysis below. No negative curvature Bianchi model admits a manifestly isotropic 
representation. 

Because of isotropy, the only non-vanishing connection one-form coefficients are nec- 
essarily of the form 

r^o^vW'^", r\^x{t)e\c<j'- (10) 

Here eabc '■— ^[abc] is the usual 3 dimensional Levi-Civita anti- symmetric symbol. 

From the definition of the curvature (|2), one can now find all the nonvanishing curvature 
2-form components: 

R\ = xdt A e\,<j- + ^2]^a ^^^^ xCe\ce'^ja' A a\ (11) 

R^o = ipdt Aa" - xi^(y^ A t'^bco" + V'Ce^cff'' A a". (12) 
Consequently, the scalar and pseudoscalar curvatures are, respectively, 

R = Q[a-^ + a-\^? ~[x~- C,f + C')], (13) 
X = Q[a-^X + 2a-^^j{x~0]. (14) 
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Because of isotropy, the only nonvanishing torsion tensor components are of the form 

T"w = u{t)5l, T\, = -2x{t)e\,, (15) 

where u and x are the scalar and pseudoscalar torsion, respectively. From the definition 
of the torsion one can find the relation between the torsion components and the gauge 
variables: 

u — a^^{a — ip), x — a^^{x — C)- (16) 

Note that a^^ip — H — u, where H = a^^d is the Hubble function. 

Regarding the material source of gravity, because of the symmetry assumptions, the 
source material energy-momentum tensor is necessarily of the fluid form with a flow vector 
along the time axis and an energy density and pressure: p, p. Although we expect the source 
spin density to play an important role in the very early universe, it is reasonable to assume, 
as we do here, that the material spin density at later times is negligible. 



4. Effective Lagrangian 

The dynamical equations for these homogeneous cosmologies could be obtained by im- 
posing the Bianchi symmetry on the general field equations found by BHN from their 
Lagrangian density. On the other hand dynamical equations can be obtained directly and 
independently (generalizing the procedure used in Ref. "SO^to include Bianchi IX, pressure 
and the new couplings) from a classical mechanics type effective Lagrangian, which in 
this case can be simply obtained by restricting the BHN Lagrangian density to the Bianchi 
symmetry (this step is where the manifestly homogeneous representation plays an essen- 
tial role). This procedure is known to successfully give the correct dynamical equations for 

T 1 

all Bianchi class A models (which includes our cases) in GR, and it is conjectured to 
also work equally as well for the PG theory. Our calculations explicity verify this property 
for the isotropic Bianchi I and IX models. Indeed the equations we obtained in this way 
are equivalent to those found (at a later date) by BHN for their FLRW models (which are 
equivalent to our isotropic Bianchi models) by restricting to FLRW symmetry their general 
dynamical PG equations. This has proved to be a useful cross check. 

Our effective Lagrangian Loff = + ^int includes the interaction Lagrangian: Li„t = 
pa^, where p — p{t) is the pressure, and the gravitational Lagrangian: 

1 3 
Lq = ■;^(aoi? + boX - 2A)a^ + — (-021*^ + Aaj,x^ + 4o-2ua;)a^ 
2k 2k 

-^{w^R^ -w^X^ + p^RX)a\ (17) 

It should be noted that the parameter restrictions 02 < 0, < 0, > 0, and + 
Aw^Wq < are, in the light of eqs. ( fTST l, (fT4l l. ( ITSI l, necessary for the least action principle, 
which requires positive quadratic-kinetic-terms. 

In the following we often take for simplicity units such that k = 1 = g. These factors 
can be easily restored in the final results by noting that in the Lagrangian they occur in 
conjunction with certain PG parameters. Hence in the final results one need merely make 
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the replacements {00,02,03,60, A, 0-2} n ^^0, 02, 03, &o. A, 0-2}, {^3, ^"6, A*3} ^ 

£I"1{W3,U;6,M3}- 

The gravitational Lagrangian has the usual form of a sum of terms homogeneous in 
"velocities" Lq = Lq + Li + L2', the associated energy function is thus 



o"^-^ — 3(ao — -a2)M'^ — 3aoH^ + 3x^(oo — 203) 



r^tp + - 




3(ao 


+6uH(ao 


W6 

~24 




w, 
24 




M3 

24 


RX 



G . . dLo . 



da 



R' - 12RUH-uY-x^ + ^ 
1^ o^ 



(18) 



The energy value ( fTSl l has the form —a^p where p can be identified as the material energy 
density. 

Making use of the formulas for the torsion and curvature components in terms 
of the gauge variables (I13I14I16] |. we now obtain the Euler-Lagrange equations 

— — — — 0, where qk = ]■;/') X: a r- The dynamical equations are the ip equa- 

dt dqk oqk 
tion: 



d OLr 



dt di> dt 

= 3(02^ — 2a2x)a^ 
M3 



3oo - —R - —X 



dip 

6ao - weR - ^X 



aip 



6bo-^R + w3X 



a(x - C), 



(19) 



which can be rearranged into the form 



^R-^X= 
2 4 



the X equation: 

d OLq 



3(2ao - 02) - - u- 6{bo + (72) - ^ R + W3X 



(20) 



dt dx dt 



dx 



-6(2030; + cr2u)o^ 



Gbo-ljR + wsX 



M3 



6ao - wqR ^X 

aip, 



a(x - C) 



(21) 
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which can be rearranged into the form 



4 2 



Q{ba + a2)-'^R + W3X 



u + 



6(ao - 2a3) - w^R - 



x: 
(22) 



[—a 3[a2U — 2a2x\) = — h 



(23) 



and the a equation: 

d OLq _ d 

dt da dt^ " "L"^- ^-i/ ' 

= 3a-' L - (I - H|r - |x) [a-R + 6(*= - [x - CI" + C')l 

+3a^{a2U — 2(72x)u — Ga^paaa; + a2u]x + 3pa^, (24) 
which can be rearranged into the form 

3 

- 3{a2U - 2a2x) = &H[a2U - 2a2x] + -{aoR + boX - 2 A) 

9 1 
+ -(-02^2 + 4a3.T2 + 4(72 Mx) - siweR^ - w^X'^ + ii^RX) 
1 8 



V 2 



12 

&0 , «'3 



24 



-R - 6{H - uf + 6x^-6 



2 ' i2^-i^n-^+i2..(i/-.)] 



+3(a2W — 2c72x)w — 6[2a3X + G2u]x + 3p 



(25) 



-3(a2U - 2(T2i;) = {aoR + 6o-'«^ - 3A) + QH[a2U - 2a2x\ 



+ -j^{—a2U + 4a3X + 4(T2Mx) 



24 



{w&R^ - w^X^ + ii^RX) 



2 
1 



/'2 

-6(il - uf + _ gS 



— {2wqR + 113X) 



+ 6box{H - u) 

-6{H - uf + 6x^ - 6^ 



+ — {2W3X - fxsR) [12x{H - u)] + 3p. 



(26) 



Since Lq is time independent, the energy function satisfies an energy conservation 
relation: 

£ = —-t-tQ = — ir-rw F — X F — a = — — Oj = 3pa a, (27) 



hence, as expected, we recover the perfect fluid relation 

d{pa^) _ da^ 
df ~^~dt' 



5a 



(28) 
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The above equations ( I20I22I26I I are 3 second order equations for the gauge potentials 
a, V', X- However they can in an alternative way be used as part of a set of 6 first order 
equations along with the Hubble relation a = aH and the following two relations, obtained 
by taking the time derivatives of the torsion (fT6t and using the curvature definitions (113114) : 



H 



-Hx 
W 



6 



- H{H -u)-{H -u) 



(29) 
(30) 



One advantage of such a reformulation is that the variables are now all observables. 

Our 6 first order dynamical equations and the energy constraint equation can now be 
put in the form 



a 
H 



aH, 

1 
6a2 



{a2R - 2a2X) - 2H^ 



a2 - 4a3 2 C 
^ — 

02 a 



3a2 



4A 
3a2 ' 



1 





u = 




X — 


--R- 


tJ-3 y 












4 



(aoi? + CT2^) - 3Hu 

3a2 

-j~i3H-2u)x, 



4a: 

02 



3 2 

X 



3a2 



3(32 + 
-6(72 



D I ^3 ^ 



/i3 



R-w-iX 



6(T2 ^ 
12^3 



Ai3 n 



W3X 



X 



(31) 



(32) 
(33) 

3a2 

(34) 
(35) 
(36) 



3(-ia2 + 253)2:2 + ^ 



- (^-6(72 



M3 



P-i 



3a2 + iy6-R+ 



-A 
1 

+ 24 

\2 „2 



(H-uY-x' + \ 



(37) 



where we have introduced certain modified pctrameteis 02, 03 and t72 with the definitions 

02 := 02 - 2ao, a3:=a3-Y, a2 := (72 + Oq. 
The last two dynamical equations (135136b can of course be resolved for R and X: 



(38) 



i? = - [(u;3a2 - P3^2)u - (2w3(T2 + 2^3a3)a;] 
a 

2a a 
(2W60-2 + -^^302)" + (4w6a3 - P3<72)x 



X = - 



-2Xu 



(4wi + Mi) 

2a 



Rx 



{W3 - W6)H3 



Xx, 



a 



(39) 



(40) 
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where 



a := 



Ms 



(41) 



Note that for the range of parameters of physical interest a > 0. 

We have cast our system into six first order equations for (3D) tensorial quantities, 
equations which are suitable for numerical evolution and comparison with observations. 
However these equations are probably not in the most suitable form for the most penetrating 
analytic analysis. So we here also present the Hamiltonian equations for our PG cosmology. 



5. Hamiltonian formulation 

From the above one can introduce the canonical conjugate momentum variables: 



Pa = = -3a^ [a2U - 2a2x] , 



P^i = 



da 
dip 



die ~ 



3ao 



36o + 



We 


^3 vl 










2 





(42) 
(43) 
(44) 



Now one can construct the effective Hamiltonian: 



^eff = PaO- + P^lp + PxX - LeS 



= a^(A-p)-6aa3(x-C) + 
Za 



2 , SaWix-0 



02 



+Pu 

+P^ 

+Px 



(f-x + c)+V' 



V-' + (x - 0' - C' + 



(60/U3 - 2aoW3)a^ 



2a 



-2V(x - C) + 
1 



(aoMs + 26oW6)a^ 



2a 



6a 



+ Pl 



W3 
6a 



Pi 





1 


, p2 


1 


1 


-6a. 


a 


^ ^ a 


602. 


a 



-■ (45) 
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From the effective Hamiltonian, we obtain the six first order Hamilton equations: 

Pa 



X 
P. 







dPa 






1 


dP^ 


a 




1 


dPx 


a 



3a2a ' 



6a 



da 



-2^(X-C)- 



At3(3a^ao - P4,) + 2we{3a%o + P^) 
6a 



— (a-X + C) + V^ 

0-2 



2a 



(A-p) 



+12a3(x-C)' + P, 



(5o// 3 - 2aoW3) 
2a 



P. 



{aofi3 + 2&0W6) 9aja{x - C) 



2a 



02 



9-H 



off 



12aa3X 



a2 a ^ 



(46) 
(47) 
(48) 



(49) 
(50) 
(51) 



dx 02 

This canonical reformulation should be of considerable interest for further studies of 
this model, since the Hamiltonian formulation is the framework for the most powerful 
known approaches for analytically studying the dynamics of a system, including such tech- 
niques as the Hamilton-Jacobi method and phase space portraits. 



6. Numerical Demonstration 

We present the results of a numerical evolution of our cosmological model. For all these 
calculations we take A = and p = 0. We need to look into the scaling features of this 
model before we can obtain the sort of evolution results we seek on a cosmological scale. In 
terms of fundamental units we can scale the variables and the parameters as k = 8ttG — 1. 
So the variables and the scaled parameters wq and W3 become dimensionless (note: the 
Newtonian limit gives oq = 1). However, as we are interested in the cosmological scale 
to see changes on the order of the age of our Universe, let us introduce a dimensionless 
constant Tq, which represents the magnitude of the Hubble time (Tq = Hq^ = 4.41504 x 
10^^ seconds). With this scaling, all the field equations are kept unchanged while the period 
T TqT. 

For our example evolution we take the parameters as 

oo = 1, 02 = —0.83, 03 — —0.35, We — —1.1, W3 — 0.091, 

(72 = 0.26 and fis = 0.21. (52) 

The behavior of the 6 equations has been observed with several sets of initial values. We 
plot this typical case in two sets of figures. In Fig.[Tl the evolution of the expansion factor 
a, the Hubble function, H, the second time derivative of the expansion factor, d, and the 
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Fig. 1 . Evolution of tlie expansion factor a, the Hubble function, H, the 2nd time derivative of the expansion 
factor, a and the mass density, p as functions of time with the parameter choice and the initial data as specified 
in the text. The gray (solid) line represents the evolution with the pseudoscalar parameters a2 and /X3 turned off. 
The red (dashed) line represents the evolution with the parameter (T2 activated. The blue (doted) line represents 
the evolution including both pseudoscalar parameters (T2 and p,^ . 



energy density, p. In Fig. |2l first line, the scalar curvature and the pseudoscalar curvature, 
R and X, second line, the torsion and the axial torsion function, u and x. 

To show the effect of the new pseudoscalar coupling parameters (T2, /is we have shown 
the evolution with these parameters turned off, with only d^, and with both pseudoscalar 
parameters activated. 



7. Concluding discussion 

We have been investigating the dynamics of the Poincare gauge theory of gravity. Recently, 
the model with two good propagating modes carrying spin 0+, 0^ (referred to as the scalar 
and pseudoscalar modes) has been extended to include pseudoscalar constants that couple 
the two different parity modes. Here we have considered the dynamics of this BHN model 
in the context of manifestly homogeneous and isotropic cosmological models. We found 
an effective Lagrangian and Hamiltonian and the associated dynamical equations. The La- 
grangian equations were rearranged into a system of 6 first order equations suitable for 
numerical evolution and a sample evolution was presented which showed the effect of the 
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Fig. 2. In the first line we compare the scalar curvature, R and the pseudoscalar curvature, X in different 
situations. In the second line we compare the torsion, u and the axial torsion, x. The first column is the evolution 
with vanishing pseudoscalar parameters, (T2 and /13, the second column, with parameter (T2, the third column, 
with both pseudoscalar parameters, a2 and ^3. 

pseudoscalar coupling constants — which provide a direct interaction between the even and 
odd parity modes. 

In these models, at late times the acceleration oscillates. It can be positive at the present 
time. It should be noted that the 0+ torsion does not directly couple to any known form of 
matter, but it does couple directly to the Hubble expansion, and thus can directly influence 
the acceleration of the universe. On the other hand, the 0^ couples directly to fundamental 
fermions; with the newly introduced pseudoscalar coupling constants it too can directly 
influence the cosmic acceleration. 

The objective of the present work was to derive certain dynamical equations for the PG 
BHN isotropic homogeneous cosmological model. These will serve as a foundation for our 
future investigations into the dynamics of this model. 

Acknowledgments 

This work was supported by the National Science Council of the R.O.C. under the grants 
NSC-98-21 12-M-008-008 and NSC-99-21 12-M-008-004 and in part by the National Cen- 
ter of Theoretical Sciences (NCTS). 

References 

1. F. W. Hehl, P. von der Heyde, G. D. Kerlick and J. M. Nester, Rev. Mod. Phys. 48 (1976) 393. 

2. F. W. Hehl in Pmc. of 6th Course of the International School of Cosmology and Gravitation 
on Spin, Torsion and Supergravity , eds. P.G. Bergmann and V. de Sabbata (New York: Plenum, 
1980), p. 5. 

3. K. Hayashi and T. Shirafuji, Prog. Theor. Phys. 64 (1980) 866, 883, 1435, 2222. 

4. E. W. Mielke, Geometrodynaniics of Gauge Fields, (Berlin: Akademie-Verlag, 1987). 



PGT with coupled even & odd parity dynamic spin-0 modes 1 3 



5. F. W. Hehl, J. D. McCrea, E. W. Mielke and Y. Neeman, Phys. Rep. 258 (1995) 1. 

6. F. Gronwald and F. W. Hehl, "On the Gauge Aspects of Gravity", in Proc. 14th Course of the 
School of Gravitation and Cosmology (Erice) , eds. P.G. Bergmann, V. de Sabbata, and H.J. 
Treder (Singapore: World Scientific, 1996), p. 148-98. 

7. M. Blagojevic, Gravitation and Gauge Symmetries, (Bristol: Institute of Physics, 2002). 

8. P. Baekler, F W. Hehl and J. M. Nester, Phys. Rev D 83 (201 1) 024001. 

9. F.-H. Ho and J. M. Nester, Poincare gauge theory with even and odd parity dynamic connection 
modes: isotropic Bianchi cosmological models, to appear in Journal of Physics: Conference 
Series, eds. L. Horwitz and M. Land (2011), arXiv: 1105 .5001 

10. H.-T. Nieh and M.-L. Yan, J. Math. Phys. 23 (1982) 373. 

11. P. Baekler and F. W. Hehl, Beyond Einstein-Cartan gravity: Quadratic torsion and curvature 
invariants with even and odd parity including all boundary terms, arXiv: 1105.3504 

12. E. Sezgin and P van Nieuwenhuizen, Phys. Rev D 21 (1980) 3269. 

13. M. Blagojevic and I. A. Nicolic, Phys. Rev. D 28 (1983) 2455. 

14. H. Chen, J. M. Nester and H.-J. Yo, Acta Phys. Pol. B 29 (1998) 961. 

15. R. Hecht, J. M. Nester and V. V. Zhytnikov, Phys. Lett. A 222 (1996) 37. 

16. H.-J. Yo and J. M. Nester, Int. J. Mod Phys. D 8 (1999) 459. 

17. H.-J. Yo and J. M. Nester, Int. J. Mod Phys. D 11 (2002) 747. 

18. H.-J. Yo and J. M. Nester, Mod Phys. Lett. A 22 (2007) 2057. 

19. K.-F Shie, J. M. Nester and H.-J. Yo, Phys. Rev D 78 (2008) 023522. 

20. H. Chen, F-H. Ho, J. M. Nester, C.-H. Wang and H.-J. Yo, /. Cosmol. Astropart. Phys. JCAPIO 
(2009) 027. 

21. A. V. Minkevich, Phys. Lett A 80 (1980) 232. 

22. A. V. Minkevich, Phys. Lett A 95 (1983) 422. 

23. A. V. Minkevich and I. M. Nemenman, Class. Quant. Grav 12 (1995) 1259. 

24. A. V. Minkevich and A. S. Garkun, Class. Quant. Grav 23 (2006) 4237. 

25. A. V. Minkevich, A. S. Garkun and V. I. Kudin, Class. Quant. Grav. 24 (2007) 5835. 

26. H. Goenner and F. MUUer-Hoissen, Class. Quant. Grav. 1 (1984) 651. 

27. X.-Z. Li, C.-B. Sun and P Xi, Phys. Rev D 79 (2009) 027301. 

28. X.-Z. Li, C.-B. Sun and P Xi, / Cosmol. Astropart Phys. JCAP04 (2009) 015. 

29. X.-C. Ao, X.-Z. Li and P Xi, Phys. Lett. B 694 (2010) 186. 

30. C.-H. Wang and Y.-H. Wu, Class. Quant. Grav 26 (2009) 045016. 

31. A. Ashetkar and J. Samuel, Class. Quant. Grav. 8 (1991) 2191. 



